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Numerical  research  is  reported  on  the  propagation  of  short  microwave  pulses  into  living,  biological  materials. 
These  materials  are  dispersive,  and  data  on  the  dielectric  constant  and  conductivity  for  these  materials  follow  a 
Debye  model.  A  Fourier-series  calculation  is  presented  that  predicts  the  occurrence  of  Brillouin  precursors  when 
the  incident  pulses  have  sufficiently  short  rise  times.  These  transients  are  attenuated  with  increasing  propagation 
distance  but  are  attenuated  more  slowly  than  the  carrier  frequency  of  the  pulse,  which  is  attenuated  exponentially 
with  distance.  An  analysis  of  the  numerical  error  resulting  from  truncation  of  the  Fourier  series  is  givsn.  Upper- 
bound  estimates  of  truncation  error  show  good  series  convergence. 


1.  INTRODUCTION 

Microwave  engineers  are  currently  developing  radiating  de¬ 
vices  that  generate  short,  high-energy  pulses.1'3  As  a  result 
it  is  a  matter  of  considerable  medical  interest  to  determine 


quency  ai.  Our  approach  is  to  treat  individual  pulses  as 
members  of  a  pulse  train,  so  that  the  analysis  is  amenable  to 
a  Fourier-series  analysis.  A  typical  trapezoidally  modulat¬ 
ed  incident  signal  (in  vacuum)  is  shown  in  Fig.  2  and  is  given 
by 


E(z,  f)  =  (l/a)(t  -  z/c)sin[oj(f  -  z/c)j 
=  sinMt  -  z/c)} 

=  |1  -  (1/a) [£  -  z/c  -  (r  -  a)]jsin[w(f  -  z/c)] 
-  0 


for  0  <  (t  —  z)/c)  <  a 
for  a  <  {t  -  z/c)  <  (r  -  a) 
for  (r  —  a)  <  (t  —  z/c)  <  t 
for  r  <  (t  -  z/c)  <  L, 


precisely  the  dynamical  evolution  of  these  pulses  as  they 
progress  through  living  tissue,  in  order  to  obtain  an  assess¬ 
ment  of  possible  health  effects  in  workers  exposed  to  these 
signals. 

In  this  paper  we  describe  some  recent  numerical  research 
on  the  dynamical  evolution  of  short-rise-time  microwave 
pulses  that  are  transmitted  into  tissue.  Since  tissue  is  com¬ 
posed  largely  of  water,  we  have  directed  considerable  atten¬ 
tion  to  the  study  of  pulse  propagation  through  this  material. 
The  frequency  dependences  of  both  the  relative  dielectric 
constant,  tr,  and  the  conductance,  c,  of  pure  water  are  shown 
in  Fig.  1.  These  curves  have  been  computed  from  Eqs.  (1.1) 
which  were  provided  to  us  by  Grant4: 

ir  =  5.5  +  72.7/[l  +  («T)2], 

a  =  1(T5  +  72.7w\T/[i  +  (air)3].  (1.1) 

These  equations  were  determined  by  fitting  experimental 
data  to  a  Debye  model  medium,  and  they  apply  for  frequen¬ 
cies  w  from  dc  to  101 1  Hz  at  a  temperature  of  25°C.  Here  e0  is 
the  permittivity  of  free  space,  T  =  8.1  X  10" 12  sec  is  the 
relaxation  time,  and  w  -  2irf,  in  which  f  is  the  frequency  in 
hertz.  For  frequencies  >10:l  Hz  we  have  used  data  that 
were  calculated  from  values  for  the  index  of  refraction  and 
the  absorption  coefficient  provided  in  Fig.  7.9  of  Ref.  5. 


2.  INCIDENT  WAVEFORMS  STUDIED 

We  have  studied  incident  pulses  that  are  formed  by  trape¬ 
zoidal  modulation  of  a  sinusoidal  field  of  fixed  angular  f re¬ 


in  which  c  is  the  speed  of  light  in  vacuum.  The  parameter  a 
is  the  rise-time  parameter  for  the  pulse  envelope  and  is 
simply  the  time  that  it  takes  for  the  modulation  to  rise  to  its 
maximum  value.  We  have  generally  assumed  that  the  maxi¬ 
mum  amplitude  of  our  incident  signals  is  1  V/m.  The  signal 
described  in  Eq.  (2.1)  is  repeated  with  a  period  L,  so  that  the 
Fourier  series  of  the  pulse  train  is 

E{z,  t)  =  (A„2  +  Bn2)1/2sin(An(t  -  z/c)  +  0„],  (2.2) 


for  which 

aLAn  =  jsin(ii)  +  X„)a  +  sin(a>  +  Xn)(r  -  a) 

-  sin(w  +  \„)rj/(a)  +  Xn)2  +  [sin(a>  -  Xn)a 

+  sin(a)  —  \„)(r  -  a)  -  sin(a>  -  Xn)r]/(a)  -  X„)2, 
aLEn  =  [1  —  cos(a)  +  X„)a  +  cos(w  +  Xn)r 

-  cos(u!  +  Xn)<7  —  a)]/(ie  +  X„)2  +  [—1  +  cos(a>  -  X„)a 

-  cos(a)  -  Xn)r  +  cos(cc  -  X„)(r  -  a)]/(u>  -  X„)2, 

where  X„  =  2im.IL  and  tan  <pn  =  AJBn. 

As  the  rise-time  parameter  a  approaches  zero,  we  find,  by 
using  L’Hopital’s  rule,  that 

limit  An  =  2co(  1  —  cos  Xn7-)L(«2  —  X„2), 

o  -0 

limit  Bn  =  —  2ai  sin  X„r /L(a/2  -  Xn:). 

a  -0 
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Fig.  1.  Relative  permittivity  or  dielectric  constant  (dashed  l’ne) 
and  conductance  (in  reciprocal  ohms  per  meter),  shown  as  functions 
of  frequency.  Logarithmic  scales  are  used. 


Fig.  2.  Example  trapezoidal  pulse.  The  trapezoidal  envelope  can 
be  modified  to  achieve  the  extremes  of  square-wave  or  triangular- 
wave  modulation. 


Therefore,  a  square-wave-modulated  sinusoidal  pulse  train 
has  the  Fourier  series 

£(z,  f)  =  V  Cn  sin[A„(f  -  z/c)  +  0„],  (2.3) 

n*=l 

in  which 

C„  =  -4w  sin(A„r/2)/L(w2  -  An2), 

<t>n  =  -A„t/2. 

The  strategy  of  our  analysis  is  to  use  these  Fourier-series 
representations  of  the  incident  pulse  train,  to  propagate 
each  Fourier  component  individually  into  the  structure  of 
interest  by  using  Maxwell’s  equations,  and  then  to  sum  the 
individual  transmitted  sinusoids  to  reconstruct  the  propa¬ 
gated  response  to  the  pulse  train.  The  pulses  in  the  incident 
pulse  train  are  assumed  to  be  spaced  widely  enough  apart  so 
that  there  is  no  mutual  interference  of  the  pulses  within  the 
target  structure. 


3  EQUATIONS  FOR  TRANSMITTED  WAVES 

The  Fourier  series  for  transmitted  pulse  trains  are  now  con¬ 
sidered  for  the  half-space  goemetry.  The  analysis  begins 
with  Maxwell’s  equations,  which  assume  the  following  form 
for  our  problem: 

V  X  E  =  -dB/dt, 

V  X  H  =  J  +  dB/dt, 

V  •  D  =  0, 

V  -  H  =  0.  (3.1) 

Here  B  =  moH  because  tissue  is  believed,  at  this  time,  to  have 
the  magnetic  properties  of  a  vacuum.  In  addition,  p  is  taken 
to  be  zero. 

In  order  to  include  the  dependence  of  a  and  t  on  the 
frequency  for  a  homogeneous,  isotropic,  locally  linear  casual 
medium  filling  the  half-space  z  >  0,  we  express  the  conduc¬ 
tion  current  J(r,  t)  in  the  constitutive  relation 

J(r,  t)  =  [  VFj(u)E(r,  t  -  iz)du,  (3.2a) 

Jo 

in  which  Wj(u)  is  the  conduction  current  impulse  response. 
Also,  the  displacement  vector  D(r.  f )  is  given  by  the  constitu¬ 
tive  relation 

D(r,  t)  *  f  WD(u)E(r,  t  -  u)du,  (3.2b) 

Jo 

in  which  Wd(u)  is  the  displacement  impulse  response. 

Applying  Maxwell’s  equations  [Eqs.  (3.1)]  with  the  above 
constitutive  relations  to  plane-wave  propagation,  in  the  z 
direction,  we  seek  solutions  of  the  form  £  =  £$(z)exp  (i\t)  and 
find,  as  is  well  known,  that  such  solutions  must  obey  the 
Helmholtz  equation, 


d2tijdz2  +  =  0, 

(3.3) 

in  which 

*2  “  A2M(jf  “  iAMotr. 

(3.4) 

Notice  that  the  permittivity  is  given  here  by 

(  =  Re(Ww)  + 

(3.5) 

and  the  conductivity  is  given  by 

a  =  Re(  Wj  -  0)  Im(tf'fj), 

(3.6) 

in  which  a  tilde  over  a  function’s  symbol  indicates  the  Fouri¬ 
er  transform  of  the  function.  These  relations  show  that  in  a 
dispersive  medium  such  as  tissue,  in  which  there  can  be 
separate  processes  underlying  the  conduction  and  polariza¬ 
tion,  (  and  a  can  each  relate  to  both  processes. 

Representing  k  as  0  —  ia,  we  obtain  the  plane-wave  phase 
factor  d  as 

=  A(Muf/2)1/2|l  +  [1  +  (<r/Xt)2],/2|l/2  (3.7) 

and  the  amplitude  attenuation  coefficient  a  as 

a  =  A(Mo(/2)1/2|-l  i-  [1  +  U/A<)2]1/2]1/2.  (3.8) 

The  half-space  geometry  considered  here  is  shown  in  Fig. 
3.  In  the  half-space  z  <  0,  the  wave  factors  «i  and  di  relate  to 
the  vacuum,  so  that  di  =  A/e  and  a;  =  0.  In  the  half-space  z 
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Fig.  3.  Diagram  of  the  half-space  geometry  studied  in  this  paper. 


>  0,  the  wave  factors  ag  and  fi n  describe  the  tissue,  so  that  «u 
and  depend  on  p  =  nn,  <  =  f(A),and<r  =  <t(M.  Application 
of  the  boundary  conditions  then  yields  the  complex  trans¬ 
mitted  field  £tr f  ,s(z,  t)  in  the  form 


(di  -  for,) 


(d;  —  in,)  +  (dn  -  ia„) 


X  exp[-f*|i?  +  f ( At  —  di|3)l,  (3.9) 


in  response  to  an  incident  wave  E ■  exp|iA(t  -  z)/c)\. 

After  taking  the  imaginary  part  of  f,rum(z.  t)  we  find  that, 
for  a  continuous  incident  wave  of  the  form  £in,(z,  t)  =  E, 
sin(A/t  -  z/c)],  the  transmitted  wave  at  depth  z  in  the  half¬ 
space  z  >  0  is  given  by 

=  2£,exp(-«„z)sin(At  -  (3uz  +  0) 
lrans  [(1  +  0nc/\)~  +  (auc/A)']1' 
in  which  the  phase  angle  C  is  specified  by 

tan  6  ~  «,,/((A/c)  +  dnl-  (3.11) 


Therefore,  if  a  pulse  train  such  as  that  given  by  Eq.  (2.2)  is 
incident  upon  the  half-space,  the  transmitted  wave  is  given 
by  the  summation 


ls(z,  t) 


-2X 


( An 2  +  Bn2)1/2exp(— an  „z)sin(Af  -  d,|  ,z  +  </>„  -i-  t)n) 


f(l  +  /3||„c/A)“  +  (a||,nc/A)‘] U1 


(3.12) 


In  this  expression  the  quantities  «n  „  and  dn,„  are  simply  an 
and  dii  calculated  at  A„.  Similarly,  0n  is  fl,  as  given  above, 
calculated  by  using  an  „  and  du.n  Tor  a  squ«re-wave-niodu- 
lated  sinusoidal  incident  wave,  the  quantity  (A„-  +  B„-)u- 
appearing  in  Eq.  (2.2)  is  replaced  by  C„. 


4.  RESULTS 

Commerfeid,”  Brillouin,'  and  Oughstun  and  Sherman*  a 
used  asymptotic  analyses  to  show  that  a  Heaviside  unit- 
step-function-modulated  sinusoidal  field  produces  funda¬ 
mental  transient  fields,  called  precursors,  in  a  medium  ex¬ 
hibiting  anomalous  dispersion.  Our  analysis  is  concerned 


with  water  and  other  similar  media  that  do  not  exhibit 
anomalous  dispersion  in  the  frequency  regime  relevant  to 
microwave  pulses.  Inspection  of  Eq.  (3.12)  shows  that  it  is 
likely  that  significant  pulse  deformation  or  blurring  will 
occur  on  propagation  into  the  medium,  but  it  is  not  evident 
‘hat  any  recognizable  transients  will  result. 

In  fact,  from  computations  of  the  transmitted  propagated 
field  with  Eq.  (3.12)  we  have  noticed  a  well-formed  transient 
field  that  appears  similar  to  the  Brillcuin  precursor  observed 
in  media  with  anomalous  dispersion.  On  the  other  hand, 
thus  far  in  our  numerical  research  we  have  observed  no 
evidence  of  the  Sommerfeld  precursor  when  we  use  dielec¬ 
tric  constant  and  conductivity  data  as  shown  in  Fig.  1.  In 
part,  this  is  because  the  Brillouin  precursor  is  a  low-frequen¬ 
cy  transient  field  and  so  i9  excited  more  readily  by  an  input 
microwave  field,  whereas  the  Sommerfeld  precursor  is  a 
high-frequency  transient  field  and  is  not  excited  so  readily 
by  such  a  low-frequency  input  field. 

Figure  4  shows  examples  of  transmitted  waves  at  increas¬ 
ing  propagation  distances  in  the  haif-space  z  >  0  when  the 
incident  wave  is  a  square-wave-modulated  1-GHz  sinusoid. 
The  transient-field  formation  is  most  pronounced  at  greater 
propagation  depths.  The  1-GHz  carrier  frequency  is  atten¬ 
uated  exponentially,  as  is  expected,  but  the  peak  transient- 
field  amplitude  decrease.”  much  more  slowly.  This  is  shown 
in  Fig.  5,  in  which  the  peak  of  the  leading  transient  is  plotted, 
along  with  the  amplitude  of  the  fundamental  frequency,  as  a 
function  of  depth  in  the  half-space.  This  behavior  is  in 
agreement  with  the  asymptotic  analysis  of  Ref.  9,  in  which  it 
is  shown  that  the  peak  amplitude  in  the  Brillouin  precursor 
decreases  only  as  z-1/2. 

In  addition,  we  have  found  that  the  transient-field  forma¬ 
tion  is  influenced  sensitively  by  the  nature  of  the  pulse 
modulation.  Specifically,  increasing  values  of  the  initial 
rise-time  parameter  a  in  the  trapezoidally  modulated  field 
results  in  decreasing  amplitudes  of  the  transient  fields 
formed,  including  a  reversal  of  the  polarity  of  the  transient- 
field  peak,  depending  on  the  size  of  a.  This  is  illustrated  in 
Fig.  6,  which  shows  plots  of  the  transmitted  signals  at  150  cm 
into  a  water  half-space  for  increasing  values  of  the  rise-time 
parameter  a. 

Our  numerical  calculations  for  10-GHz  square-wave-mod¬ 
ulated  signals  produce  pulse  waveform  patterns  that  are 
strictly  analogous  to  those  shown  in  Fig.  4  for  the  1-GHz 
pulses  with  the  exception  of  a  simple  scaling  by  the  penetra¬ 
tion  distance.  For  example.  Fig.  7  shows  a  plot  of  the  tran¬ 
sient  peak  amplitude  and  the  amplitude  of  the  fundamental 
frequency  versus  the  propagation  distance,  and  this  is  seen 
to  be  similar  to  that  depicted  in  Fig.  5  when  the  distance  axis 
is  scaled  linearly. 

The  numerical  results  displayed  in  Figs.  4-7  were  ob¬ 
tained  by  direct  numerical  calculation,  using  Eq.  (3.12).  For 
each  frequency  A„  of  the  Fourier  series  we  used  a  value  of 
<(A„)  and  a  value  of  <x(A„)  taken  from  a  spline  fit  of  the  data  in 
Fig.  1.  In  the  case  of  the  1-GHz  pulse,  r  was  10  nsec  in 
duration  and  L  was  20  nsec  in  duration.  In  the  case  of  the 
10-GHz  pulse,  t  was  1  nsec  and  L  was  2  nsec.  We  have  found 
in  all  our  calculations  thus  far  that,  when  L  g  2r,  individual 
pulses  in  the  pulse  train  do  not  interact  and  therefore, 
when  the  pulses  are  adequately  spaced,  one  obtains  simply 
the  individual  pulse  response. 
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Fig.  4.  Tiansmitted  wave  at  depths  of  (a)  5  cm,  (b)  75  cm,  and  (c) 
150  cm  in  a  water  half-space  when  the  incident  wave  is  a  square- 
wave-modulated,  1-V/m,  1-GHz  sinusoidal  field  orthogonally  inci¬ 
dent  upon  the  air-half -space  interface. 


Fig.  5.  Peak  of  the  leading  transient  (dashed  curve)  and  amplitude 
of  the  fundamental  frequency  (solid-curve),  as  functions  of  the 
propagation  depth  in  the  half-space  z  >  0  for  the  1-GHz,  square- 
wave-modulated  signal. 


We  have  calculated  the  transmitted  wave  by  using  Eq. 
(3.12)  and  summing  over  the  coefficients  from  n  -  1  to  and 
including  n  —  512.  Truncation  error,  caused  by  failure  to 
sum  terms  associated  with  n>  512  can  be  evaluated  directly 
from  Eq.  (3.12).  For  example,  for  a  square-wave-modulated 
incident  wave,  the  neglected  signal  E513  (z,  t)  is  given  by 


^513(2,  t) 

=  2  V  Cn  eXP(~“H."Z)sin(Xt  ~  A'!."2  +  ^  +  0n) 

n^3  K1  +  0II.aC/X>2  +  («I1,-.C/X)2]1/2 

Since  aii,n  and  du,n  are  both  always  greater  than  or  equal  to  0, 
and  since  exp(— «n,nz)  and  sin(u)  are  always  less  than  or 
equal  to  1, 


4u) 


8u> 


R, ,n(z,  t)  s  2  V  — - - —  <  • 

A  W  ~  «2)  WW-«*) 


8w  r  d  n 

L  J513  (X„2  -  w2) 


(4.2) 


It  then  follows  that  the  neglected  signal  satisfies  the  inequal¬ 
ity 


< 


8u> 

^(X5l32  —  <‘>2) 


X513  ~  <*> 
X513  +  w 


(4.3) 


This  last  expression  demonstrates  that  the  series  for  the 
transmitted  wave  converges  for  all  z  g  0.  For  the  1-GHz 
square-wave-modulated  pulse  studied  in  this  report. 
^5i3(z»  t)  <  0.0248  V/m.  Improved  estimates  of  R 5i3(z,  t)  for 
z  >  0  can  be  calculated  by  carrying  the  attenuation  term 
exp(— aii,„z)  in  the  estimate.  By  doing  this,  one  finds  a 
progressive  decrease  in  the  truncation  error  with  increasing 
propagation  distance.  As  empirical  evidence,  when  we  use 
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(a) 


Fig.  6.  Transmitted  wave  at  a  depth  of  150  cm  in  a  water  half-space 
for  an  incident  wave  that  is  a  trapezoidally  modulated,  1-V/m,  1- 
<1Hz  sinusoid.  The  rise-time  parameter  a,  the  time  to  reach  the 
maximum  value  of  the  trapezoidal  modulation,  is  set  at  (a)  1/(2 f) 
sec,  (b)  3/(4 f)  sec,  and  (c)  3/(2 ft  sec. 


Fig.  7.  Peak  of  the  leading  transient  (dashed  curve)  and  the  ampli¬ 
tude  of  the  fundamental  frequency  (solid  curve),  as  functions  of  the 
propagation  depth  in  the  half-space  z  >  0  for  the  10-GHz,  square- 
wave-modulated  signal. 


1000  or  more  coefficients  to  calculate  transients,  we  general¬ 
ly  note  amplitude  changes,  compared  with  the  500-coeffi¬ 
cient  calculation,  in  only  the  fourth  or  fifth  decimal  place. 
By  comparison,  although  a  similar  approach  using  the  fast- 
Fourier-transform  algorithm  significantly  increases  the 
speed  of  calculation  and  seems  to  work  quite  well,  we  have 
not  yet  been  able  to  develup  an  estimate  of  the  error  of  the 
calculated  transmitted  wave  because  of  the  finite  sampling 
of  the  incident  signal. 

Experiments  have  been  performed  at  Physics  Internation¬ 
al  Company  in  San  Leandro,  California,  to  test  the  above- 
described  theory  of  pulse  propagation.  Pulse  propagation 
in  a  coaxial  test  cell  and  in  an  open-field  tissue-equivalent 
tesi  structure  has  followed  the  numerical  predictions  closely. 
This  research  will  be  reported  fully  elsewhere. 

5.  CONCLUSION 

We  conclude  that  a  short-rise-time  microwave  pulse  can 
create  a  disturbance  in  tissue  at  a  greater  depth  than  would 
be  expected  if  one  were  to  use  th.  depth  of  penetration  of  the 
basic  carrier  frequency  «.  This  finding  should  be  of  interest 
to  those  who  are  studying  or  otherwise  evaluating  human 
exposures  to  electromagnetic  devices  that  produce  such 
pulses.  In  addition,  from  our  calculations  of  square-wave- 
modulated  sinusoidal  signals,  we  have  observed  a  transient 
field  in  water  that  is  similar  to  the  second  or  Brillouin  pre¬ 
cursors  that  were  described  asymptotically  by  Brillouin7  and 
by  Oughstun  and  Sherman,6’9  Finally,  we  have  observed 
that,  for  trapezoidal  modulation  of  the  input  field,  this  tran¬ 
sient  field  decreases  in  amplitude  as  the  time  to  reach  the 
maximum  of  the  trapezoidal  modulation  is  increased,  and  an 
inversion  of  the  transient  polarity  is  produced. 
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